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ABSTRACT 

Im;; 'irities enter a snowpack as a result of fallout or scavenging by falling snow crystals. Al- 
bedo and flux extinction coefficient of soot contaminates snowcovers are studied using a two-stream 
approximation o!' the radiative transfer equation, The ef'ect of soot is calculated by two methods: 
independent scattering by ice grains and impurities, and average refractive index for ice grains. Both 
methods predict a qualitatively similar effect of soot: the albedo is decreased and the extinction co- 
efficient is increased compared to that for pure sn<iv/ in the visible region; the infrared properties 
are largely unaffected. Quantitatively, however, the effect of soot is more pronounced in the average 
refractive index method. We find that soot contamination provides qualitative explanation for sev- 
eral snow observations. 
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ALBEDO AND FLUX EXTINCTION COEFFICIENT OF IMPURE 
SNOW FOR DIFFUSE SHORTWAVE RADIATION 


Albedo and flux extinction coefficient, which are basic shortwave radiative properties of 
snow, are known to be highly variable. For deep snowcovers observed albedo ranges from 65 to 
96 percent and extinction coefficient from 0.07 to 0,17cm'i in the blue-green (0.4Sjum) region 
of the solar fiux (Thomas, 1963: Kondratiev et al,, 1964; Grenfell and Maykut, 1977; Mellor, 
1977), Recent radiative transfer studies (Choudhury, 1981; Wiscombe and Warren, 1980) could 
not explain such variability for pure snow. Since attempts are being made to correlate changes 
in albedo and microwave brightness temperature with snowpack physical characteristics (McGinnis, 
1975; Foster et al„ 1980), it is important to understand this variability. 

Natural snowcovers are expected to contain some impurities. Soot, released into the atmos* 
phere due to fossil fuel burning, can enter snowpacks by fallout or as a result of scavenging by 
falling snow_^ crystals. The absorption coefficient of soot is more than five and eight orders of 
magnitude higlier than pure ice in, respectively, microwave and visible regions (Bartky, 1968). A 
small concentration of soot can therefore affect visible and microwave observations. Althougli 
soot is not the only impurity that can be present in a snowcover, this report gives quantitative re- 
sults for soot only, because of its high absorptivity. 

Multiple scattering within the snowpack is studied by a two-stream approximation, for which 
Choudhury and Chang (1979) and Grenfell and Maykut (1977) have given equations for albedo 
and flux extinction coefficient. The scattering parameters (cross-sections and the asymmetry 
factor for phase function) needed for the calculation of these radiative properties depend upon 
size and refractive index of scattering particles. For impure snow, consisting of particles (ice 
grains and impurities) with differing refractive indices, the scattering parameters can be calculated 
by two methods (Bartky, 1968; Lindberg, 1975; Gillespie et al., 1978; Warren and Wiscombe, 
1980): (1) calculate separately the sc.attering parameters for ice grains and impurities and aver- 



age by corresponding number fractions, and (2) use fractional volumes of ice and impurity to cal- 
culate an average refractive index for all particles. These methods, referred below as isolated par- 
ticle and average refractive index models, may represent physical situations wlien tlie impurities 
are unattached or attached to the ice grains. Certainly, these two situations cail co-exist within 
a snowcover. A recent study by Warren and Wiscombe (1980) is baseO on the isolated particle 
model, and they used a different radiative transfer solution. 

The scattering parameters are calculated using geometrical optics formulae (Choudliury and 
Chang, 1981) for ice grains and Mie scattering theory (Hansen and Travis, 1974) for soot parti- 
cles. We assumed all ice grains are of equal size, and used a modified-gamma size distribution 
for soot particles (Bartky, 1968). The refractive indices are from Hobbs (1974) for ice and from 
Dalzell and Sarofim (1969) for propane soot. The soot parameters were chosen sligliUy different 
by Warren and Wiscombe (1980). We note that the characteristics of the impurity will depend 
largely upon its source. 

Albedn of pure and soot contaminated snow are shown in Figure 1. Soot in amounts of a 
few parts per million by weight (ppmw) can reduce visible albedo by several percent without 
affecting near infrared albedo. It is also seen that although isolated particle and average refrac- 
tive index models predict similar effects of soot, this effect is more pronounced in the average 
refractive index model, It is difficult to clarify In physical terms the cause of this difference 
between the moiels since they use different methods for calculating the scattering parameters. 
However, if these models are representing the physical situations discussed previously, this differ- 
ence indicates that the impsmities imbedded or attached to the ice grains are more effective in 
changing the albedo than the isolated impurities. For a given weiglit fraction of soot, the cor- 
responding volume fraction increases as the ice grain size increases. Thus for same weight frac- 
tion, soot will affect the visible albedo of course grained snow more than fine grained snow. If 
new snow is contaminated (e,g., due to scavenging) then with aging the visible albedo will decrease 


rapidly due not only to increasing grain size but also increasing volume fraction of impurity. 

These consideraTiuns are pertinent in remote sensing when snow depth is to be correlated with 
visible albedo. 

Comparison with observations (Grenfell, 1981) is shown in Figure 2. The ice grain radius 
of 0.2mm, used in this comparison, is slightly larger than the observed 0.1 to 0.15 mm. Al- 
though by varying the amount of soot the visible albedo can be matched using thr, observed grain 
radius, the near infrared albedo will be predicted a few percert higher. This discrepancy witli 
respect to grain size indicates some deficiencies in the two-stream radiative transfer model of 
snow. One should, however, note that near infrared albedos are affeijted primarily by the ice 
grains in the surftice centimeter layer. If these grains are slightiy larger that fire observed values 
(due to noniiniforra metamorphism) then this discrepancy will not arise. The grain profile In- 
forraation is not given In Grenfell (1981) to further analyze this possibility, 

Bohren and Barkstrom (1974) showed that tire albedo is directly proportional to the square 
root of grain size. We calculated integrated albedos for varied grain size and soot content, imd 
found such a relationship. For the average refractive mde-K model, the intercept (C| ) and slope 
(C 2 ) of this relationslup, A * Cj - Cir^^ (where r is grain radius in mm), are given in Table I. This 
parameterized equation may be helpful in climatology and hydrology for energy balance calculations. 

In Figure 2 a fairly good agreement between observed and calculated albedos is shown. A 
crucial model validation, however, is in matcliing both albedos and flux extinction coefficients, 
when they ha/e been observed concurrently. The reason being they provide complementary radia- 
tive transfer information dependent on the same parameters. One should, however, note that 
the extinction coefficient depends more strongly on grain siie and absorption coefficients of ice 
and impurhy than does albedo. An attempt to match such observations by Grenfell and Maykut 
(1977) is shown in Figure 3. An ice grain radius of 0.1mm seems more appropriate for matching 
tire albedos because for this radius the predicted albedos agree better witli observations beyond 
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O.fJSMm - a region not affecteti by soot. With this grain radius if we adjust the soot covdent so 
as to match the visible albedo, we find that the predicted extinction coefncients disagree by more 
than a factor of two. it is also seen that matching of extinction coefficients was not very .suc- 
cessful, The weakness of the model and/or inaccuracies of input parameters (both soot and ice) 
to represent these observations are clear. Warren and Wiscombe (1980) matched these albedos 
but did not compare with extinction coefficients. In matching the albedo observations of Kuhn 
and Siogas (1978), Warren and Wiscombe used l.Sppmw of soot with ice grain radius O.i mm. 

For these soot and ice parameters we find that calculated extinction coefficients are an order of 
magnitude higher than the observations.. In fact tlie extinction coefficients observed by Kuhn 
and Siogas are what may be expected for pure snow (Choudhury, 198!). It is possible that these 
discrepances tme parUy due to observational error (M. Kuhn, private communication) and partly 
due to assumed impurity characteristics, but untii explained, the model predictions should be 


treated cautiously. Our results, however, substantiate tile importance of impurity contamination 
(Dunkle and Bevans, 1956: Warren and Wiscombe, 1980) in understanding the radiative proper'* 
ties of natural snowpacks, 


A part of the work was done while B. J. Choudhury was at Computer Science Corporation 
under Task Assignment 125 of NASA contract NAS-24350. 
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Figure 1 . Effect of 1 .0 ppmw of soot on albedo of snow with grain radius 0.1 mm. (A) pure snow, 
tB) isoKated particle model, (O average refractive inde,N: model. For the same amount of soot note 
the difference between (B) and tC). 
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Figure 2. Comparison of calculated albedo with observations (Grenfell. 1981). (A) pure snow, (B) 
isolated particle model with l.Oppmw of soot, (C) average refractive index model with O.Sppmw of 
soot. Note the difference in the amount of soot lor (B) and (C), 
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Figure Captions 


Figure I, Effect of LOppmw of soot on albedo of snow with grain radius 0.1 mm. (A) pure 
snow, (B) isolated particle model, (C) average refractive index model. For the same 
amount of soot note the difference between (B) and (C). 

Figure 2. Comparison of calculated albedo with observations (Grenfell, 1981), (A) pure snow, 
(B) isolated particle model with l.Oppmw of soot, (C) average refractive index model 
with 0.5 ppmw of soot. Note the difference in the amount of soot for (B) and (C). 

Figure 3, Comparison of flux extinction coefficient and albedo with obsers'ations (Grenfell and 
Maykut, 1977). Results from isolated particle and average refractive index models are 
shown side-by-side. Note the difference in soot content for these results and the 
sensitivity of extinction coefficient with respect to ice and soot parameters in the 
visible region. 
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